Soliton solutions of 3D Gross-Pitaevskii equation by a potential control method 



Renato Fedele, 1 ']*] Bengt Eliasson , 2 ■ 3 ' f \ Fernando Haas, 3, [^] Padma 
Kant Shukla, 3 ' 4 f]Dusan Jovanovic, 5 Q and Sergio De Nicola 6 '0 

1 Dipartimento di Scienze Fisiche, Universita Federico II and INFN, 
Complesso Universitario di M. S. Angelo, Via Cintia, 80126 Napoli, Italy 
2 Department of Physics, Umea University, SE-90 187 Uraed, Sweden 
3 Institut fur Theoretische Physik IV and Centre for Nonlinear Physics, 
Fakultat fur Physik und Astronomie, Ruhr-Universitdt Bochum, D-44 780 Bochum, Germany 
^Scottish Universities Physics Alliance (SUPA), Department of Physics, 
University of Strathclyde, Glasgow G4 ONG, United Kingdom 
5 Institute of Physics, P. O. Box 57, 11001 Belgrade, Serbia 
6 Instituto di Cibernetica "E Caianiello" del CNR Comprensorio 'A. Olivetti", 
Via Campi Flegrei 34, 1-80078 Pozzuoli (Na), Italy 
(Received 8 July 2009) 

We present a class of three-dimensional solitary waves solutions of the Gross-Pitaevskii (GP) equa- 
tion, which governs the dynamics of Bose-Einstein condensates (BECs). By imposing an external 
controlling potential, a desired time-dependent shape of the localized BEC excitation is obtained. 
The stability of some obtained localized solutions is checked by solving the time-dependent GP 
equation numerically with analytic solutions as initial conditions. The analytic solutions can be 
used to design external potentials to control the localized BECs in experiment. 

PACS numbers: 02.30.Yy, 03.75.Lm, 67.85.Hj 



I. INTRODUCTION 



Solitons are nonlinear localized wave packets sustained by the balance between wave dispersion and medium non- 
linearity. Solitons propagate over large distances without changing their shape [TJ [51 [3] . Solitary waves or solitons 
have been observed in several areas of physics including fluids, plasmas, optics, biology, and condensed matters (e.g. 
Bose-Einstein condensates). Many types of solitons have been studied, starting with classical examples found in 
integrable models, such as the Korteweg-de Vries, sine-Gordon, Toda-lattice, and nonlinear Schrodinger equations, 
and their non-integrable extensions. Solitons are robust against collisions due to the integrability of the underlying 
equations. 

Recent observations of matter-wave solitons [H [5j [6j [7j |8j |9j [101 E] have been among the most groundbreaking 
achievements in the burgeoning fields of Bose-Einstein condensation (BEC) of dilute atomic gases. In the latter, 
bosonic atoms below a certain temperature suddenly develop in the lowest quantum mechanical state. The balance 
between the spatial dispersion of matter waves and repulsive or attractive atomic interactions in Bose-Einstein con- 
densates (BECs) ensures the existence of dark j4j|5l|6j|7] or bright [8, 9, [TO] [TTj solitons, respectively. A dark (bright) 
matter wave soliton is a localized BEC having a minimum (maximum) condensate density at the center. However, 
if the atomic condensate is embedded into a periodic potential created by standing light waves, i.e. optical lattice 
|12j . there exists possibility of reversing the matter wave group dispersion sign (e.g. from positive group dispersion 
to negative group dispersion), and of the possible observation of the bright matter- wave soliton in the BECs with 
repulsive inter-atomic interaction. The concept of the dispersion control by periodic potentials is also well known in 
solid state physics |13j and a very active topic of research in nonlinear optics |14j . The dynamics and stability of the 
matter wave solitons in the BECs is governed by the nonlinear Schrodinger equation [T], known as the Gross-Pitaevskii 
(GP) equation jT5j [T6] in the context of the BECs. For BECs with positive (repulsive) interparticle interactions, dark 
solitons with locally depleted density have been studied theoretically |17j and have been observed in many experi- 
ments jH |5J O [71 IHJ 121 EH EC] • Numerical and theoretical studies revealed the dynamics and stability the BECs in 
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a magneto-trap both for repulsive and attractive interactions in a limited parameter range in spherically symmetric 
trap [18] . The stabilization and controlling of the BECs in asymmetric traps have been investigated by considering 
the time-dependent solutions of the GP equation [5D]. Stable condensates with a limit number of atoms of 7 Li with an 
attractive interaction have been observed in a magnetically trapped gas |19j . The formation of matter- wave solitons 
[U [5] and trains of solitons [SJ H] have been observed in BECs of 7 Li atoms that are confined in a quasi-one-dimensional 
trap and magnetically tuned from repulsive to attractive interactions. The solitons are predicted to cither collapse 
or explode, depending on the parameters of the BECs and on the confining or repulsive potential |21j . The collective 
collapse |22l I23j and explosion |23j of BECs with attractive interactions have been observed in experiments. Several 
theoretical and experimental studies of coherent matter waves are contained in Ref. [23]. Wang et al. [25] have 
presented the analytical dark and bright solitons of the one-dimensional GP equation with a confining potential. The 
generation of matter wave dark and bright solitons in a prescribed external potential for confining the BECs have 
been investigated with a periodically varying nonlinear coefficient [26, 27J. A tight transverse trap with a gradually 
varying local frequancy along the longitudinal direction induces an effective potential for one-dimensional solutions in 
a self-attracted BEC [28]. The propagation of a dark soliton in a quasi-lD BEC in the presence of a random potential 
has been studied by Bilas and Pavloff [29 . Furthermore, there is a recent theoretical study of exact one-dimensional 
solitary wave solutions in a radially confining potential |30j . The dynamics of the one-dimensional bright matter 
wave soliton in a lattice potential has been studied by Poletti et al. Two-dimensional dark solitons to the 

nonlinear Schrodinger equation are numerically created by two processes to show their robustness (e.g. stable against 
the head-on collision). Stcllmcr et al. [32] present experimental data exhibiting the head-on collision of dark solitons 
generated in an elongated BEC. Experiments do not report discernible interaction among solitons, demonstrating the 
fundamental theoretical concepts of solitons as quasiparticles. 

In this paper, we theoretically study possibility of controlling the time-dependent dynamics of the BECs in three- 
dimensions with a carefully spatially shaped, time-dependent controlling potential. This idea is theoretically supported 
by our recent mathematical investigations [33] and on the improvement of the recently-proposed 'controlling potential 
method' [34]. From the experimental point of view, this idea could be realized by techniques involving lithographically 
designed circuit patterns that provide electromagnetic guides and microtraps for ultracold systems of atoms in BEC 
experiments |35j . and by optically induced "exotic" potentials |36J. The stability of the obtained analytic solutions 
are investigated numerically by direct integration of the time-dependent GP equation. 



II. THEORY 

The dynamics of BECs is in a non-uniform potential is governed by the three-dimensional Gross-Pitaevskii equation 
[I3CL6] 

ih— = -—V 2 * + gN\*\ 2 * + 14xt(r, , (1) 
at 2m 

where m is the atomic mass, and V ext is the external potential, g = 4irh 2 a/m where a is the short-range scattering 
length, which can be either positive or negative, giving rise to either repulsive or attractive interactions, and N is the 
number of atoms in the BECs. Typical parameters values where solitons have been observed in BECs of 7 Li atoms 
are N = 10 4 -10 5 at a temperature of 1-10 /xK and a magnetic field ~ 400-600 G, leading to a small scattering length 
of a w — 0.2nm [U [9]. In Eq. (Ill, W is normalized such that J l^l 2 d 3 r = 1. Equation ^ can be cast into the 
dimcnsionless form 

^^ = ~V 2 V> + 5|V'lV + ^Wr,^)^, (2) 

where the time t is normalized by to = (47r) 2 ma 2 iV 2 /Ti, space r by ro = 47r|a|iV, the external potential £7 ex t by 
h 2 /(47r) 2 ma 2 N 2 , and ip = fQ 2 so that J \ip\ 2 d 3 r = 1 in the normalized spatial variables. With this normalization, 
we have g = +1 for a > and g = — 1 for a < 0. For typical experimental values N = 6000, \a\ — 0.2 nm and 
m = 1.17 x 10 -26 kg ( 7 Li) [5J, we would have ro = 16 (im and to = 0.028 s. Our goal is to design the external potential 
to obtain a desired time-dependent shape of the solution. 

We concentrate on a sub-class of solutions, in which the GP equation can be separated into one linear, two- 
dimensional equation, and one nonlinear, one-dimensional equation. In doing so, we first make a decomposition of 
the external potential according to |33] 



U ext (r,t) = U ± (r ± ,t) + U z (r ± ,z,t), 



(3) 
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where rj_ is the position vector perpendicular to the z direction, and make the ansatz that the solution can be 
separated as 



where the normalization conditions j \ ipx\ 2 d rx — 1 and f \ip z \ 2 dz = 1 are imposed. 
Inserting Eqs. Q and Q into Eq. ([2]), and reordering the terms, we have 



(4) 



i^ + lv 2 x iPx-Ux(rx,t)iPx 



-tpj 



(5) 



where we have denoted = d 2 jdx 2 + d 2 /dy 2 . By requiring that tpx satisfies the linear Schrodinger equation 



r + lVx^-Ux(rx,t)i;x=0, 



dt ' 2 

we obtain the one-dimensional GP equation 

+ YW^ ~ ^lUlV* - U z (r ± , z,W, = 0. 
For the linear Eq. ([6]) for ij)x (r_L > t), we assume a parabolic potential well in the form 



(6) 



(7) 



Ux = ^K x {t> 



-K v {t)y 2 



(8) 



which yields solutions of Eq. (16]) in terms of functions of the form tpxmn(x,y,t) = ipxm(x,t)ip yn (y,t) where ip a p 
(where a = x, y and (3 = 0, 1, .77) are Hermite-Gauss functions in the form 



, _ exp[-a 2 /4q^)] 

M , ; " [2^2 (t)2 2/3 (/ 3 ! )2]l/4 11 /5 



a 

V2a a (t) 



exp 



« 2 1" *<Pa/j(*) 



0) 



(normalized so that \ip a /s(a,t)\ 2 da = 1) where H^^) are Hermite polynomials of order /3. The first few Hermite 
polynomials are listed in Table [T] in the Appendix A. Here the phases are given in terms of a a as 



7a (t) 



1 C?<7 Q 



and 4>a/3(t) = (2/3+ l)</> Q0 (i), where 



d<t>aa{t) _ 1 

~ 4o£(f)' 



We have that a a is related to K a by the Pinney equation [39 



(10) 



(11) 



(12) 



Solving for K a (t) in (12) and inserting the result into (18j) we have the potential well 



l / l l\ 

X 2 V ct x di* ^ Aa*) 



1 d 2 a y 1 \ 



(13) 



in terms of (Ta;(£) and (J y (t). Hence, we have the possibility to arbitrarily choosing the time-dependent widths a x (t) 
and <J y (i) (and the indices m and n) of our solution, and obtain the potential U perp necessary to sustain that solution. 
Now, the solution of Eq. ^ for ip z must be such that the total ip(r,t) in Q solves the original GP equation (|2]). A 
special solution is 



ip z (z,t) 



gS m 5„ 



Aa x (t)a v {t) 



1/2 



sech 



2<J X {t)(Jy{t)' 



exp 



-g(t)z 2 + i6 (t) 



(14) 
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where phase functions are given by 



and 



1 



9 = 



d[a x (t)a v (t)} 



a x {t)a v (t) 



dt 



d6a 
dt 



1 



g5 m 8 n 



_a x {t)a y {t) 

The external potential U z to sustain this solution is given by 



U z (rj_, z,t) 



\1>±\- 



a x (t)a y (t) 



(15) 



(16) 



(17) 



where 



K(t) = 



<J X {t)<Jy{t) d^ 



[<T X (t)<Ty(t)} 



(18) 



The details of the derivation of i/) z and U z are given in the Appendix A. The perpendicular solutions ip±_ in ([9| and 
parallel solution ip z in (14 1 can now be used to construct the total, three-dimensional solution ip in Q of the GP 
equation and where the external controlling potential U ext in ^ is the sum of U± in (13 1 and U z in ( fl7| . 



III. ANALYTICAL AND NUMERICAL EXAMPLES 




FIG. 1: The soliton widths a x (solid line) and a y (dashed line) as functions of time, during one oscillation period, for the 
equilibrium widths <jq x = ao y = 0.075, the oscillation amplitudes ao x = ao y = 0.25, and the frequency lo — 300. 



We will here exemplify our analytic results with exact time-dependent solutions of the GP equation. As mentioned 
above, we are free to choose arbitrary time-dependencies of the widths a x (t) and u y (t). In Eq. (|9j) we also need the 
Hermite polynomials H m and in Eqs. (pL4|) , (|16| and (|17|) we need the constants <5 m , which are given in Table [I] m 



Appendix A, for different excited states m. As an example, we choose a x (t) and cr y (t) to be periodic in time, of the 
form 



and 



& x {t) 



<Ty(t) 



&0i 



1 + ao x sin(ujt) 



<Jq v 



1 + a 0y cos(ujt) 



(19) 



(20) 



where we investigate cases of relatively large amplitude oscillations with ao x = ao y = 0.25, and we set a$ x — <7q v — 
0.075 and u> — 300. In dimensional units with r = 16 /im and t = 0.028 s, this corresponds to a typical soliton width 
of - 0.075 x 16^m « 1.5 ^m and the frequency 300/(2tt x 0.028) Hz = 1.7kHz. The choice of the functions ([l9l 
and (20 1 describes periodically in time pulsating solutions widths, illustrated in Fig. [lj Small values of a x and a y 
correspond to spatially localized solitary waves, while large values correspond to wider solitary waves. The amplitude 
of the solitary wave also varies so as to keep the total number of condensates constant, J |-0|e? 3 r = 1. In Figs, pj-pl 
we have plotted the solutions together with the corresponding controlling potentials in the xy-plane at z = 0, for 
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FIG. 2: The condensate density \ip\ 2 (left) and the external potential U ex t (right) in the xy plane at z = 0, at times t = 0, 
t — 0.25 x 2n/u>, t = 0.5 x 2-k/ui, and t = 0.75 x 2tt/lj (top to bottom rows), for the ground state (ra, n) = (0, 0). 




FIG. 3: The condensate density |i/>| 2 (left) and the external potential U ex t (right) in the xy plane at z = 0, at times t = 0, 
t = 0.25 x 2n/u>, t = 0.5 x 2-k/cj, and £ = 0.75 x 2-k/uj (top to bottom rows), for the excited state (ra, n) — (1, 0). 

the ground state (m, n) — (0,0) and the excited states (m,n) = (1,0) and (m,n) — (1,1), where (m,n) refers to the 
orders of the Hermite polynomials in the perpendicular solution ipj_ mn (x,y 1 t). The widths are varying such that at 
time t = 0, the solitons are localized and large amplitude while at later times their amplitudes decrease and their 
widths increase, first in the y direction and then in the x direction. We note that the external potential is sometimes 
large amplitude and confining, and sometimes small amplitude and non-confining. 

For the analytic results to be observable in experiments, it is necessary that they are stable. To assess the stability 
of the solutions, we have therefore solved the time-dependent GP equation (2) in three dimensions with the analytic 
solution as initial condition at time t = 0. We used a box length L x = L y = 2 in the x and y dimensions and L z = 3 in 
the z direction, with periodic boundary conditions. The spatial derivatives were approximated with a pseudospectral 
method, and the time stepping was performed with the standard 4th-order Runge-Kutta method. The space was 
resolved with 64 grid points in the x and y directions and with 200 grid points in the z direction, and the timestep 
was taken to be At = tt/300 000. 

In order to seed any instability in the system, we added random perturbations in phase to the initial condition 
of the order 0.01 rad. We then simulated the system for 5 oscillation periods and measured the maximum density 
\ip\num max obtained in the numerical solution as well as the maximum relative error in density fluctuations, e — 
|(|V'| 2 — \4 , \num)\max/\'4'\ m axi as a function of time, and plotted the results in Fig. [5] for the ground state (m, n) — (0, 0) 
and the excited states (m,n) = (1,0) and (m, n) — (1,1). We see that the numerical solution follows almost exactly 
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FIG. 4: The condensate density |i/>| 2 (left) and the external potential U ex t (right) in the xy plane at z = 0, at times t = 0, 
t — 0.25 x 2n/u>, t = 0.5 x 2-k/ui, and t = 0.75 x 2tt/lj (top to bottom rows), for the excited state (m, n) — (1, 1). 



Density Relative error 
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FIG. 5: Numerical simulation results of the GP equation for the amplitude ao x — ao y — 0.25 and (m,n) = (0,0), (1,0) and 
(1,1) (top to bottom panels), showing the maximum density (left column) and the relative deviation of the numerical solution 
from the exact analytic solution (right column). 



the analytic solution, without increasing substantially throughout the simulation. The ground state (m, n) = (0,0) 
shows a relative error less than 1 % throughout the simulation while the excited state (m, n) = (1, 1) shows a somewhat 
larger error but still less than 10 %. Hence the controlled BEC seems to be stable enough to be observed in experiment. 



IV. SUMMARY AND CONCLUSIONS 



In summary we have presented a class of three-dimensional solitary waves solutions of the Gross-Pitaevskii (GP) 
equation, which governs the dynamics of Bose-Einstein condensates. This has been done on the basis of our recent 
mathematical investigations [33 and improving the formulation of the recently-proposed controlling potential method 
[34] . By imposing an external controlling potential, a desired time-dependent shape of the localized BECs is obtained. 
The stability of the exact solutions were checked with direct simulations of the time-dependent, three-dimensional GP 
equation. Our simulations show that the localized condensates are stable with respect to perturbed initial conditions. 
We propose that our findings could be tested experimentally by techniques involving lithographically designed circuit 
patterns that provide electromagnetic guides and microtraps for ultracold systems of atoms in BEC experiments [35j , 
and by optically induced "exotic" potentials |36j . Furthermore, numerical simulations |40j reveal that soliton emission 
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from a BEC can be controlled by a shallow optical dipole trap. Here, the emission of matter wave bursts is triggered 
by spatial variation of the scattering length along the trapping axis. The motion of the ID dark soliton can also be 
controlled by means of periodic potentials in optical lattices [41 . Finally, it should be noted that repulsive BECs 
confined by an optical lattice and a parabolic magnetic trap can appear in the form of vortices |42] as well. Computer 
simulations |43j reveal the condensation of a finite temperature Bose gas in the form of a single vortex. Interactions 
of solitary waves and vortex rings in a cylindrically controlled BECs exhibit robustness during head-on collisions |44j . 



APPENDIX A: DETAILS IN THE DERIVATION OF THE PARALLEL SOLUTION ip z AND THE 

POTENTIAL V z 

We here present a special solution of Eq. Q for ijj z , such that the total Tp{v,t) = ip±(r±,t)ip z (z,t) solves the 
original GP equation S. Multiplying Eq. ftjy by IV'jJ 2 and integrating over r_L space, we obtain 



i^+l^T- «o(t)l^| 2 ^ - V(z,t)il>, = 0, 



where V(z, t) = J U z (r±, z, t)\ijj±\ 2 d 2 r± and 



Qo(t) = g / dxdy = 



gS m S„ 



where the numerical factor 



a x (t)a y (t) ' 

1 f°° 
^—^1 exp(-2^)|H„«)r <J{ 



is evaluated for given values of m. The first few Hermite polynomials and values of 6 m arc listed in Table ||J 



(Al) 



(A2) 



(A3) 



m 


Hm(0 


5m 





1 


1/(V5F) 


1 


2£ 


3/(8Vi) 


2 


4£ 2 - 2 


41/(128^) 



TABLE I: The Hermite polynomial // m (^) and the value of 5 m for m — 0, 1, and 2. 

Equations ^ and ( |A1[ ) should give the same solution for ip z , and this imposes restrictions on the external potential. 
Subtracting Eq. ([7]) from Eq. ( Al I we obtain the compatibility relation 



U g {t ±1 z,t) = -fl|^x| 2 |V z | 2 + <Zo(*)|V>*| 2 + V. 



(A4) 



We will consider the special case where the external potential for the one-dimensional GP equation ( Al I has the form 
V(z,t) = (l/2)K(t)z 2 so that 



. dj) z 1 d 2 tp z 
1 dt 2 dz 2 



q (t)\i> z \ 2 ip z -*K(t)z 2 il> z = 0. 



(A5) 



Using the Madelung fluid ansatz [37] ip z = y/ p(z, t) exp[i9(z, t)] , we obtain after some manipulations [cf. Eq. (11) of 
Ref. [38]] 



and 



where 



dv dp 
P dt +V i + 2 



dp d , . 

d-t + dz {pv) = °> 



co(t) 



dv 

di 



dz 



dp 

dz 



dU nTT dp 

p lTz +2U 



1 d 3 p 
dz J ' 4 dz 3 



= 0, 



d6 
dz' 



(A6) 



(A7) 



(A8) 
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In Eq. (A6), co(t) is an arbitrary function of t to be determined below and we have denoted U = qo(t)\tp z \ 2 + V(z, t) 
qo(t)p + V(z, t). In order to find solitary wave solutions of Eqs. (A6l and (A7l, we now assume that 



v(z,t) = g(t)z, 



(A9) 



where g(t) (not to be confused with the coupling parameter of the GP equation) is a function to be determined, and 
we have from Eq. ( A8 ) that 



(z,t) = -g(t)z 2 + 9 (t), 



(A10) 



where 9o(i) is an arbitrary function of t. Even if a time-dependent phase is included in the transverse solutions ip± 
below, it is necessary to keep 0q, since the parallel solutions are required to satisfy a separate equation. Using Eq. 
(A9), together with the relations dv/dt — g'(t)z and J (dv/dt)dz = (1/2) g' (t)z 2 (where the primes denote derivatives) 
in Eqs. (A6l and (A7l, we have 



g'(t)zp + g(t)z^ + 2 



^QoPw- - [ Kzp+Kz 
dz \ dz 



dp 

dz 



2 dp\ t l d 3 p 
Adz 3 



(All) 



= 



and 



m +9(t)^ + g(t) P = 0, 



(A12) 



respectively. Multiplying Eq. ( A12 1 by g(t)z and subtracting the result from Eq. (All), and reordering the terms, 
we have 



dp 



z[g'(t) + g 2 (t) + K(t)] [z^f + p) + 2c Q (t) -f - 3q oPl f + - 



t) 



dp 



dp 1 d 3 p 



dz 



dz Adz 3 



If g(t) obeys the Riccati equation 



then Eq. (Al3l simplifies to 



g'(t) + g 2 (t) + K(t) = 0, 



^ i \ dp o dp 1 d 3 p 
dz dz 4 dz 6 



We now introduce the change of variables £ — G(t)z + R(t) and r = t, or 

z = [t-R(T)]/G(r) 

and 

t=T, 



which can be introduced into Eqs. ( A13 1 and (A12) to obtain 



dp , 1 G 3 d3 P 



and 



respectively. Choosing 



G'(r) 
G 



g(r) 



dp j. n'M 9p 
M + R{T) dt 



i?'(r)^+.gp+^ = 0, 



dp 
dr 



(A13) 
(A14) 

(A15) 

(A16) 
(A17) 

(A18) 
(A19) 



G'(t)/G(t)+ 9 (t) = 0, 



(A20) 
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and R ^constant (R = without loss of generality), Eq. (A 19 1 becomes 

dp 



Ot 



(A21) 



We now look for a solution in the form p(£,r) = A(t)F(£). With this ansatz, Eqs. ( A18 1 and (A21) can be written 
as 



and 



2co(T)f(0 - 3 go (r)A(r)F(r)F'(0 + \g\t)F"'{0 = 0, 



A'(T) + g(r)A(r) = 0, 



(A22) 



(A23) 



respectively. For consistency, Eq. (A22) for F(t) should be reduced to an equ ation with the co efficients independent 
of t, i.e. cq(t) and qo(T)A(r) must be proportional to G 2 (r). Equations (A20) and (A23l also imply that A is 



proportional to G. Integrating Eq. (A20l as 



G(T) = G exp(- g(s)d i 



where we without loss of generality can choose Go = 1, we then have 

A(t)=A q G(t), 



and 



cq(t) = G G 2 (r), 



<Z„(t) = QoG(t), 



so that Eq. (A22l takes the form (after eliminating the common exponential factor) 

2C F'(0 - 3A Q F'(QF(0 + = 0. 



(A24) 

(A25) 
(A26) 

(A27) 
(A28) 



and AqQq < in the form 



Equation (A28l, which is the time-independent Korteweg-dc Vrics equation, admits solitary wave solutions for Go < 

m " (I) < (A29 > 



where A = 1/^2\C \. It follows that 



p(z,t) = ^G(i)sech 2 



Gft)z 
A 



(A30) 



with Ao 7^ 0, and where we have chosen Qq = — 2/A so that go(i) = ~(2/A)G to ensure that f_ \ip z \ 2 dz 

gS m S n 



J^pdz = 1. From (A2| we then have 



Since G = 1 at t — 0, we find the width 



so that 



A 



Q — 

A a x (t)a y (t)' 
2a x (0)a y (0) 



G 



<T X (P)<Ty(0) 

<J X (t)(Jy(t) 



(A31) 



(A32) 



(A33) 
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and 



9 = 



1 

a x (t)a y (t) 



d[a x (t)a y (t)) 
dt 



Finally we have 



g8 m 5 n 



Aa x {t)a y {t) 



1/2 



sech 



gS m S n 
2a x (t)a y (t)' 



exp 



-g(t)z 2 +i9 (t) 



(A34) 



(A35) 



where g is given by (A34l. To obtain an expression for 6>p( t) in terms of a x and a y , we use Eq. (A35l to calculate 
U(z,t), p(z,t) and v(z,t), substitute the result into Eq. (A6l, taking into account that g = —G/G. This gives 
d8 /dt — — c , which together with (A26l, (A32l and (A33l yields the result 



d6 Q 



1 



gS m Sn 



dt 8 \_a x {t)a y {t) 



Using (A2| into (A4|, we obtain the external potential 

U z (r ± ,z,t) = -g \ib ± \ 2 - 



a x (t)a y (t) 

where K is found from the Riccati equation (A14l and Eq. ( A34| as 

1 d 2 



\^ Z \ 2 + -K(t)z 2 , 



K(t) 



<7x(t)cr y (t) dt 2 



(T X (t)(Jy(t)} 



(A36) 



(A37) 



(A38) 
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